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Abstract 

Using the dinuclear system concept we present calculations of production 

cross sections for the heaviest nuclei. The obtained results are in a good 

agreement with the experimental data. The experimentally observed rapid 

fall-off of the cross sections of the cold fusion with increasing charge number 

Z of the compound nucleus is explained. Optimal experimental conditions 

for the synthesis of the superheavy nuclei are suggested. 
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I. INTRODUCTION 



In order to reach superheavy elements and the island of stability at Z = 114 and N = 
178 — 184, two heavy nuclei must fuse. At the GSI (Darmstadt) the elements with Z =110, 
111 and 112 were recently synthesized in cold fusion reactions [|l|. The heaviest isotope of 
the element with Z = 110 was produced in the FLNR, JINR (Dubna) |fj]. The nucleus with 
Z = 110 was also produced in LBL (Berkeley) ||. The next important step is the synthesis 
of the elements with Z =113 and 114 by using both Pb-based and actinide-based reactions 

§• 

The aim of investigations of the fusion mechanism is the choice of the optimal experimen- 
tal conditions for the formation of the superheavy elements. Producing the elements from 
Z = 104 to Z = 112 in the cold fusion reactions the experimentalists observed the rapid 
fall-off of the evaporation residue cross section (about four orders of magnitude) with in- 
creasing Z of the compound nucleus [|I],|],||,|7] . The measured cross section of the production 
of the element with Z = 112 is only a few pbarn. 

For the cold fusion produced elements with Z = 104 — 1 12 it was found that the cross 
sections are maximal at energies below the Bass barrier with excitation energies of the 
compound nuclei of 9-15 MeV ]3]j§|. The macroscopic models in literature HHTOfl do not 



reproduce the excitation function for these reactions. The optimal excitation energy of 
the compound nucleus is much smaller than the energy predicted in models taking a large 
extra-extra push energy into account |10 1 . In this paper we explain these effects by a fusion 



model [|ll],|r^] based on the dinuclear system (DNS) concept |L3|] and give estimations of the 
complete fusion probability and optimal excitation energy for the production of compound 
nuclei between Z = 104 and Z = 114. 

In the DNS model ||TT|--|i"3|j the fusion process is assumed as a transfer of nucleons from 
the light nucleus to the heavy one. The DNS evolves as a diffusion process in the mass 
asymmetry degree of freedom 77 = (A\ — A<i)j A to the compound nucleus {A\ and A2 are the 
mass numbers of the nuclei and A = Ai + A 2 ). Evolving to the compound nucleus the DNS 
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should overcome the inner fusion barrier Bf us in the mass asymmetry degree of freedom. 
The top of this barrier (the Businaro-Gallone point at r\ — t\bg) coincides with the maximum 
of the DNS potential energy as a function of rj. We assume that complete fusion occurs after 
the DNS overcomes this inner barrier. In the DNS-concept the value of Bj us represents a 
hindrance for complete fusion of the initial DNS with \t]i\ < \tibg\- Besides the motion in rj 
a diffusion process in the variable of the relative distance R between the DNS nuclei occurs 
leading to a decay of the DNS which we denote as quasifission. For quasifission, the DNS 
should overcome the potential barrier B q f which coincides with the depth of the pocket in 
the nucleus-nucleus potential. The energy required to overcome the fusion and quasifission 
barriers is contained in the excitation energy of the DNS. 

In the DNS concept the chosen potential energy surface does not much change during 
the time needed for the fusion through the Businaro-Gallone maximum in the mass asym- 
metry coordinate. This maximum exists in both our and adiabatic considerations. If the 
system reaches this maximum in the diffusion process, fusion inevitably occurs much easier 
in asymmetrical systems (r] > 0.75). Since in fusion reactions a gradual transition between a 
frozen density (multinucleon transfer reactions) and an adiabatic density (fission) happens, 
the complete quantitative understanding can be reached only with a cumbersome time- 
dependent calculation with a time-dependent multidimensional potential surface. In order 
to estimate the fusion cross sections in the reactions leading to the superheavy elements in 



a simple manner, we use the DNS-model JllHl3| because this model was successful in the 
description of the fusion of heavy nuclei. 

II. MODEL 
A. Evaporation residue cross section 

In accordance with the DNS-concept the evaporation residue cross section can be written 

as 
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J max 



&ER(E c , m .) — ^ °c(-Ec.m., J)PcN(E c .m., J)W sur (E c , m _, J). 
J=0 



(1) 



The value of J ma x depends on E cm _ and is smaller than Js f =o at which the fission barrier in 
the compound nucleus vanishes. Since the superheavy nuclei do not exist with large angular 
momentum, we can further use the following factorization [EJ: 



ffER[Ec.m.) — &c{E c . m )PcN{E c .ni)W sur {E cm ). 



(2) 



The capture cross section <j c defines the transition of the colliding nuclei over the Coulomb 
barrier and the formation of the DNS when the kinetic energy E cm _ is transformed into the 



excitation energy of the DNS. In the calculation of a c = cr c (E c m , J) ~ tcX 2 (J„ 

j=o 



hl) 2 T 

(A is the reduced de Broglie wavelength and T the transmission probability through 
the Coulomb barrier) for reactions leading to superheavy nuclei, the values J ma x = 107i and 
T = 0.5 are used for E rrn near the Coulomb barrier. The values of a r obtained with this 



expression are in agreement with the ones obtained in a microscopical calculation |T4j based 
on the model ||15|| . The probability of complete fusion Pcn depends on the competition 
between the complete fusion and quasifission processes after the capture stage in the DNS. 
The surviving probability W sur estimates the competition between fission and neutron evap- 
oration in the excited compound nucleus. The competition between the complete fusion and 
quasifission is not considered in the existing macroscopical models p|-|i~0"|| . 

Dissipative large-amplitude collective nuclear motions, which occur in fission, quasifis- 
sion, fusion and heavy- ion reactions, can be analyzed within the transport theory |Ll| . If the 
initial DNS (rj = rji) is at the local minimum of the potential energy in R and r/, we can use 
a two-dimensional Kramers-type expression (quasistationary solution of the Fokker-Planck 



equation) |T7j for the rates of fusion (k = if) and quasifission (k = R) through the fusion 
(B v = BJ US ) and quasifission {Br = B q f) barriers: 
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With these rates the fusion probability can be calculated 
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\Kr \Kr\Kr _ _ 

p = \ \ *R T v -Tr 

where (3 = e — 1 ~ 1.72. The first term in @) yields the contributions of the quasistationary 



rates. The second term is related to the transient times T/-. It was shown in |12j that 
we can neglect this term for <C 1/Af rr (k = R,rj) or tr ~ r v . This is true for all 
reactions under consideration. The detailed discussion of Eq. ([|) is given in Ref. ||12|| . The 
application of the Kramers-type expression |18|] to relatively small barriers (Br/® > 0.5) 



was demonstrated in [19]. The local thermodynamic temperature G is calculated with the 



expression 6 = ^fE*/a, where a = A/12 MeV" 1 and E* is the DNS excitation energy. In 
Eq. (||), the frequencies oof k ' (k,k' = R,r]) of inverted harmonic oscillators approximate 
the potential in the variables R and rj around the tops of the barriers Bk>, and Uk are 
the frequencies of the harmonic oscillators approximating the potential of the initial DNS. 
Since the oscillator approximation of the potential energy surface is good for the reactions 
considered, we neglect the nondiagonal components of the curvature tensors in (^). The 
friction coefficients are simply approximated by 7^/ = Y^yjh, ||20|| . The quantity T denotes 
an average double width of the single-particle states. The calculation of the mass parameters 
Hrr and [i m is given in [^l|] where it is demonstrated that one can neglect the nondiagonal 



mass coefficient [Ir^ in the DNS for \r/\ < \t]bg\- As shown in |T2|], the friction coefficients Jrr 
and 7 W obtained with T = 2 MeV have the same order of magnitude as the ones calculated 
within other approaches. 

The motion of the DNS to smaller rj leads also to quasifission because the quasifission 
barrier in R decreases quickly with rj due to the increasing Coulomb repulsion. As in 
Refs. |TT|-|r3|, we can use the quasifission barrier B q f for the initial DNS in (^) with a good 
accuracy. In the considered reactions the initial DNS is in the local minimum in rj and B q f 
is smaller than the barrier for the motion to smaller values of r\. 
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B. Potential energy 



The values of frequencies in ([|) are easily calculated with the potential energy of the 



DNS 111,12 



U(R, V , J) = B 1 + B 2 + V(R,r), J) - [B 12 + V rat (J)\. (5) 

Here, B±, B 2 , and B\ 2 are the realistic binding energies of the fragments and the compound 
nucleus ||22|| , respectively. The shell effects are included in these binding energies. The 
isotopic composition of the nuclei forming the DNS is chosen with the condition of a N/Z- 
equilibrium in the system. The value of U(R, T], J) is normalized to the energy of the rotating 
compound nucleus by Bi 2 + V^, ot . The nucleus-nucleus potential V(R, rj, J) in fl5|) is calculated 
as described in [|TT1 , |i"2f . Due to the large moments of inertia of the massive DNS considered 
and due to the restricted set of angular momenta (J < lOti), we can neglect the dependence 
of U(R,r],J) on J: U{R,r},J) « U{R,r}), V(R,r],J) « V{R,rj). The calculated driving 
potential U(R m ,r]) = U{rj) as a function of r\ for the reaction 54 Cr+ 208 Pb is presented in 
Fig. p] for J = 0. For a given t], the value R m coincides with the position of the minimum of 
the potential pocket in V(R,rj). 

Deformation effects are taken into account in the calculation of the potential energy 
surface [P~2"|,|2"5|] . For the heavy nuclei in the DNS, which are deformed in the ground state, 



the parameters of deformation are taken from Ref. |24]]. The light nuclei of the DNS are 
assumed to be deformed only if the energies of their 2 + states are smaller than 1.5 MeV. 
As known from experiments on subbarrier fusion of lighter nuclei, these states are easily 
populated. For the collision energies considered, the relative orientation of the nuclei in the 
DNS follows the minimum of the potential energy during the evolution in 7]. We find that 
the values of Pcn calculated with the deformation of both nuclei in the DNS are practically 



the same as the ones calculated previously in 1 12,23] where a deformation only in the heavy 



nucleus of the DNS was taken into account. However, a deformation of both nuclei in the 
DNS yields better agreement with the experimental excitation energies of the compound 
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nucleus for the In reactions. Taking the deformation of the nuclei in the DNS as a function 
of rj, we choose the way for the DNS evolution on the potential energy surface calculated in 
the DNS concept. 



III. RESULTS AND DISCUSSIONS 

In all considered cold fusion Pb-based reactions the dependence of the potential energy 
of the DNS on mass asymmetry is similar to the one presented in Fig.l. The evolution of the 
DNS in n to the compound nucleus is accompanied by overcoming the inner fusion barrier 
with the height Bj us . Due to the deformation of the DNS nuclei, which should be taken 
into account for the In reactions, the value of B*j- us decreases (the value of Pqn increases) 
as compared to the case of spherical nuclei. The values of Bj us and B q f are changed from 
4.8 and 4.0 MeV, respectively, for the 50 Ti+ 208 Pb^ 258 104 reaction till 9.0 and 1.0 MeV, 
respectively, for the 70 Zn+ 208 Pb^ 278 112 reaction. 

It is important to determine the excitation energy E^ N corresponding to the maximum 
of the excitation function in the In fusion reaction. In order to minimize the loss because 
of the fission of the excited compound nucleus, the excitation energy should be kept as 
small as possible. With the value of the inner fusion barrier, the optimal excitation energy 
is calculated as Eq N = U(R m ,7]i) + B*j- us (U(R m) rji) is the energy of the initial DNS). 
Therefore, the optimal kinetic energy is E c . m . = E^ N — Q which excesses the entrance 
barrier by the value AE = E^ N — B q f. Note that all considered collisions occur above 
the calculated entrance barrier (AE > 0). For smaller and larger excitation energies, the 
evaporation residue cross sections in the In fusion reaction decrease due to the decrease 
of the values of Pcn and W sur in @. The calculated optimal values of E^ N (see Fig. 2) 
are in good agreement with the experimental data on the In fusion reactions used in the 
production of the heaviest elements with Z =102-112 [|]]. The macroscopical models PJTO]] 
overestimate the minimal values of E^ N in the fusion reactions. The model |]TU[ predicts 
Eqpj between 50 and 300 MeV. With the model [[J the values of Eq N are estimated to be 
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about 40-50 MeV. The use of the Bass potential overestimates the experimental value of 
E* CN by 5-7 MeV §. 

The calculated values of Pcn for the In Pb-based reactions are presented in Fig.3. 
These values are in agreement with the ones extracted from experimental data [|I],|]||[7| . 
The decrease of Pcn i n (@) D Y about four orders of magnitude with Z increasing from 
104 to 112 explains the observed rapid fall-off of the evaporation residue cross sections. 
The factors a c and W sur do not strongly change with Z for the cases considered. The 
fusion probability strongly decreases with decreasing mass asymmetry of the initial DNS 
(increasing Z for the Pb-based reactions) because the inner fusion barrier Bj us increases 
and the quasifission barrier B q f decreases (the Coulomb repulsion increases). For example, 
in the 76 Ge+ 208 Pb— >- 284 114 reaction the estimated value of oer is near the limit of present 
measurements. The probability to obtain the nucleus with Z = 116 in the 82 Se+ 208 Pb 
reaction is smaller than this limit. From our analysis it follows that the fusion of symmetric 
combinations = 0) for the synthesis of the heaviest elements yields smaller cross sections. 

Using the data in Fig.3, we can explain the smaller yield of the nucleus with Z = 110 
in the reaction with 62 Ni than the one with 64 Ni. The fusion probability in the reactions 
with 66 > 68 Zn is larger than the one with 70 Zn. However, W sur in the reaction with 70 Zn can 
be larger than W sur in the reactions with other Zn isotopes because of the smaller neutron 
separation energy in 278 112. It could be that the increase of Pcn is compensated by a 
decreasing W sur in the reactions with the lighter isotopes. In addition, to obtain the same 
values of a c for the reactions with 70 Zn and 68 Zn, the excitation energy in the reaction with 
68 Zn should be larger by 2 MeV than the one with 70 Zn [p~JI]. This means that for reactions 



with lighter isotopes the optimal excitation energy obtained in the static calculation and 
presented in Fig. 2 could be within 2 MeV smaller than the realistic values. However, this 
deviation is within the present experimental accuracy. We note that in the reactions used 
for the production of the heaviest elements all factors in (Q) are equally important. 

In order to calculate the evaporation residue cross sections in the In Pb-based reactions, 
we use values of surviving probabilities W sur (E CN ) ps r n /r/ (the values given refer to an 
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angular momentum of zero) which are few times larger than the ones estimated in ||, but 
smaller than the values from the analysis of 4n reactions |[25|. In accordance with the 



experimental data and shell-model calculations |26j the value of Y n /Y f increases slightly for 
Z = 108 because of the shell closures in the vicinity of iV = 162. Since for larger Z the 
neutron separation energies and fission barriers are almost the same |2B| for the compound 
nuclei in Table 1, we took the same value of T n /T f for these nuclei. The value of T n /T f is 
sensitive to shell effects and excitation energies and has to be studied in further details. The 
calculated values of &er (Table 1) are in a good agreement with the known experimental 
data [l],[|J|[7|]. One can see that at fixed W sur and small change of a c the value of o\ n 
decreases by two order of magnitude from the nucleus with Z = 108 to the nucleus with 
Z = 113 due to the decrease of Pcn- Therefore, in the reaction 74 Zn+ 208 Pb— > 282 114 we 
expect a value of 0\ n smaller than 0.1 pbarn. 

In recent experiments the present limit of the heaviest element production in the cold 
fusion has been reached. More asymmetric combinations of the colliding nuclei than in the 
Pb-based reactions (the initial DNS is near or behind the top of the inner fusion barrier in 
mass asymmetry) can be used to extend the production of superheavy elements. According 
to our model one should take targets heavier than Pb. For the actinide-based reactions 
with the projectiles like 48 Ca, 34 ' 36 S, the fusion probability is much larger than in Pb-based 
reactions. This effect can compensate the increase of the fission of the compound nucleus 
due to a higher excitation energy which corresponds to the 3n — 4n channels. Our calculated 
cross section for the 4n reaction 48 Ca+ 244 Pu— > 292 114 is about 1 pbarn. 



IV. SUMMARY 

In conclusion, the calculated results for the (HI, In) reactions used in the production of 
the heaviest elements are in agreement with the experimental data. The calculations for all 
reactions were performed with one set of parameters and with the same assumptions. The 
main factor which prohibits the complete fusion of heavy nuclei is the quasifission. Without 
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regarding the quasifission, the explanation of the experiments on the fusion of heavy nuclei 
is not possible. Our model is useful in calculating the optimal excitation energy and the 
combinations of the colliding nuclei. It should be applied for the further analysis of the 
experimental data. 
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TABLES 



TABLE I. Calculated [th.] and experimental [exp.] evaporation residue cross sections for 
several In Pb-based reactions. The values of E c . m . lead to the values of Eq N presented in Fig. 2. 
The values of a c , Pcn and W sur used in the calculations are explained in the text. 
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FIGURES 

FIG. 1. Dependence of the potential energy of the DNS on r/ for the 54 Cr+ 208 Pb reaction 
= 0.59). The calculated results with and without deformation of the DNS nuclei are presented 
by solid and dashed lines, respectively. The result with deformation is presented for rj- values which 
are of interest to calculate Bj us . 

FIG. 2. Optimal excitation energies of the compound nucleus for the In Pb-based reactions. 
The nuclei with even and odd Z are produced with 208 Pb and 209 Bi targets, respectively. The 
experimental data are shown by solid diamonds. The projectiles are indicated. The experimental 
point for Z = 112 is shown for 70 Zn as a projectile. The calculated results are depicted by open 
circles for different projectiles. The values of Eq N obtained with the Bass potential are presented 
by the dashed line [4]. The solid line is drawn to guide the eye. 

FIG. 3. Calculated fusion probability Pcn for different In Pb-based reactions. The projectiles 
are indicated. The excitation energies in the calculations are taken the same as in Fig. 2. 
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Fig. 1. G.G.Adamian et al. 
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Fig. 3. G.G.Adamian et al. 
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